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Context. In the upper layers of the solar atmosphere the temperature increases sharply. We studied possibility of the transfer of neutrals 
motion energy into the electrostatic waves.Electrostatic waves could damp in the upper layers of the solar atmosphere and their energy 
could be transformed into the thermal energy of the solar atmosphere plasma. 

Aims. We have used the two fluid approximation when studying plasma dynamics in the low altitudes of the solar atmosphere. In 
order to study evolution of disturbances of high amplitudes the parametric resonance technique is used. 

Methods. The dispersion relation for the electrostatic waves excited due tot he motion of neutrals is derived. The frequencies of 
electromagnetic waves which could be excited due to existence of the acoustic wave are found. The increment of excited electrostatic 
waves are determined. 

Results. The motion of the neutrals in the lower solar atmosphere, where ionization rate is low, could excite electrostatic waves. In the 
upper solar atmosphere the ionization rate increases and motion of the neutrals could not support electrostatic waves and these waves 
could damp due to the collision of the charged particles. The energy of the damping waves could be transformed into the thermal 
energy of the plasma in the upper atmosphere. 
Conclusions. 
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1. Introduction 

In the solar atmosphere the temperature raises sharply from pho- 
tospheric 6000K to few 10 6 K in the corona. The problem of heat- 
ing solar atmosphere is a long standing problem in solar Physics. 
The detailed physical mechanism of coronal heating is not yet 
well understood, fundamental problems have to be solved in or- 
der to find proper theoretical models explaining observational 
data. Fundamental questions related with the coronal heating in- 
clude identification of the source of thermal energy, proposition 
of the energy release mechanism and comparison of theoretically 
predicted values of the physical characteristics with observations 
(Klimchuk, 2006). Investigation of the coronal heating needs the 
study of the highly disparate spatial scales, in particular the study 
of MHD wave dispassion, MHD reconnection, quasi-periodic 
oscillations and small-scale bursts occurring in corona is nec- 
essary (Klimchuk, 2006, Walsh & Irelan, 2003).In addition, the 
heating of the chromosphere requires much more energy, and 
this related process is equally difficult to understand. 

Based on the UV and X-ray observations of the sun, Parker 
suggested (1988) that the X-ray solar corona is created by 
the dissipation at the many tangential discontinuities arising 
spontaneously in the bipolar fields of the solar active regions 
as consequence of random continuous motion of the foot- 
points of the field in the photospheric convection. Parker named 
these elementary energy release vents as nanoflares. Number 
of studies were dedicated to finding observational evidences 
of nanoflares (Cargill, 1994, Sakamoto et al., 2008, Vekstein, 
2009). Prominent researchers found Parkers's idea very promis- 



ing and investigated nanoflares contribution into energy budget 
of the solar corona (Klimchuk, 2006). 

Parker's model does not address following questions: where 
the excess magnetic energy comes form and what is its original 
source? According to Rappazzo et al., (2008) slow photospheric 
motions induce a Poynting flux which saturates by driving an 
anisotropic turbulent cascade dominated by magnetic energy. In 
this case magnetic field lines are barely tangential and current 
sheets are continuously formed and dissipated (current sheets 
correspond to tangential discontinuities Proposed by Parker). 

Another mechanism was suggested by Krasnoselskikh et al., 
(2010), that electric waves could be generated by motion of neu- 
tral gas. A huge energy reservoir exists in the form of turbulent 
motion of neutral gas at end beneath the photosphere. It is widely 
accepted that this energy could be partially transformed into the 
excess magnetic energy in the chromosphere and corona. 

The mechanism of transformation of neutral gas motion en- 
ergy into the electromagnetic disturbances was investigated by 
Aburjania & Machabeli (1998). In this paper authors derived 
dispersion equations for coupled acoustic and electromagnetic 
modes. The authors considered studied parametric instability. 
This technique enables to study nonlinear processes, i.e. pump- 
ing energy from the acoustic waves with high amplitudes. In 
Aburjania & Machabeli (1998) authors implied discussed mech- 
anism to explain observed data and predict characteristic ex- 
pected atmospheric wave preprocesses before the earthquakes. 

In order to study excitation of the electromagnetic waves in 
solar atmosphere we use data, describing parameters of solar at- 
mosphere given by Vernaza et al. (1981). We consider the case 
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where excited velocities, electric field and wave vector and are 
parallel to external magnetic field. 

2. Main consideration 

In the low layers of the solar atmosphere (up to 1000 km, see 
Vernaza et al. 1981), the concentration of neutral is much higher 
than concentration of the charged component of plasma. At this 
altitudes of the solar atmosphere the motion of neutrals could 
influence dynamics of electrons and ions, therefore the energy 
of motion of neutrals could be transformed into the electromag- 
netic energy, in particular energy of the acoustic waves could be 
pumped into the energy of the electrostatic waves. The inten- 
sive acoustic wave draws charged particles into the motion due 
to friction with neutrals. 

Here we investigate how the energy of acoustic waves with 
high amplitude, could be transformed into the energy of electro- 
static waves. When considering the disturbances of high ampli- 
tudes, we could not neglect nonlinear terms, more over these 
terms play a key role in the process. The nonlinear pumping 
of the electromagnetic wave energy into the plasma, due to 
the parametric instability was studied by Silin (1973), in or- 
der to investigate neutral wave energy pumping into the plasma 
the parametric instability mechanism was studied by Aburjania 
Machabeli (1998), in this work we also apply parametric insta- 
bility technique. 

If we write and sum the equations of motion for electrons 
and ions, 

we will get: 

dv e dvi 
meUe ~dt +miH '~dt = 
-V (p e + pd - niiUiVin (vi - v n ) + niiriiVie (v e - t/,) 

+m e n e v ei (Vi - v e ) + m e n e v en (v n -v e ). (1) 

Here is assumed that charged particle move along the mag- 
netic field lines ant the direction of electric field and wave vector 
coincides with the direction of the magnetic field, i.e we consider 
electrostatic waves. 

In the low altitude (beyond 1000km) of the solar atmosphere 
the relations given below by Eq. 2a-2d are satisfied: 



(2a) 
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(2b) 



(2c) 



If the electric field changes with frequency much less than 
the cyclotron frequency we could assume that the charged par- 
ticle moves in a constant electric field during the Larmour ro- 
tation period. In this case the velocity of the charged particle 
is: v az = q a —t (here indices a shows the sort of particle and t 
is time). In order to define the interval between collisions of the 
charged particle kind a with neutrals we use the law of statistical 
distribution (see Aburjania & Machabeli 1998), i.e the probabil- 
ity that collision of particle a with neutrals happens between t 
and t + dt is: f an = v an e~ v ""'dt. Multiplication of v az = q a ^f-t by 



fan — Va 



(2d) 



here: 

yS = ^f- . Eq. 2d helps us to solve us three component hy- 
drodynamic equations system. 

After taking into account relations given by Eq. 2a-2d and 
comparing similar terms in Eq. (1), we obtain: 



ffij/tj dvi/dt riiOjp 



m e n e dv e /dt n e a> 



- » 1, 



pe 



mniVin 



» 1, 



m e n e v en V III 
m e n e v ei (v e - vi) 

Ve » V t . 



» 1, 



(3a) 
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(3c) 



(3d) 



After taking into account Eq. (3a-3d) and neglecting electron 
inertia, we could rewrite Eq. 1 as follows: 

% + (Vv i )v i = 
at 



1 



■ V (p e + Pi) ~ V in (Vi - V„) + V ie (V e - Vi) . 



(4) 



(5a) 



niitii 

If we represent ion velocity as: 

vt -v„ = -Vi, 
we express the velocity of neutrals as follows: 

v„ = V cos(Qt + x i'), (5b) 

Eq. 5b describes oscillatory motion of the neutrals, this type 
of motion accompanies the excitation and propagation of waves 
interacting with neutrals. In Eq. 5b *P is casual, initial phase of 
neutrals oscillations. Because of the collisions between particles 
the phase of the oscillations changes casually, the final expres- 
sions of the variables obtained after solving hydrodynamic equa- 
tions should be averaged by x ¥. 

After substitution of (5a-5b) in Eq. (4) we get 



dVi 

dt 



y O sin(Qf + ¥) + ikViVo cos(Qr + *P) 



V (p e + Pi ) - Vi (v in + pv ie ) - v ie /3V Q cos(Of + T), (6a) 

niiiii 

and 

d V 

— --V Q [Q sin(Qf + Y) + v ei cos(Of + x ¥)]+ikV i V () cos(Qr+T0 = 
at 



- — V (p e + pi) - Vi (V in + V ei ) ■ 

niitii 



(6b) 



It is convenient to use notation: 

Vt = U exp [ia sin(Qf + ¥)] . (7) 
Substitution of Eq. (7) into Eq. (6b) leads to equation: 



'dt and integration by time leads to the following q\j 



relation between velocities of electron and ions: 



- — V sin(Of + *¥) + v ei cos(Qf + T)] exp [-ia sin(Qf + ¥)] 
dt 
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— V (p e + Pi ) exp [-ia sin(Of + ¥)] - U (v in + v ei ) . (8) 

niiiii 

Let us write equation of continuity: 

d(m + n e ) 



dt 



+ V(«o + n e )v e = 0, 



(9a) 



unperturbed value of concentration-no does not depend on 
time, therefore -Jp = 0. It is convenient to introduce notation 
n e = N e (t)exp[-ibsm(Q.t + W)], here a = ^,b = 0^,0 = 



ypf. For N e we get equation: 

dN e 



dt 



ikj3U(t)n Q exp [-ib sin(Qf + *¥) + ia sin(fif + T)] , (9b) 



solution of Eq. 9b gives 



n e = ikno/3 exp[/ft sin(Qf + TOJx 

U(co') exp [-i(cj' + tjQ.)] 



X £./„(&) expO^F) J 



For pressure gradients we have: 



co' + rjQ, 



1 



m,no 
1 



-v Pi 



k z V: 



mm 



-Vpe = 



2 T /2 
Ti_ 

CO 

mm 



dco'. (9c) 

(10a) 
(10ft) 



here V Ti = -fe 

After substitution of expressions (9c), (10a), (10b) into equa- 
tion (8) and averaging by *F (it is casual phase of neutrals oscil- 
lations) we could derive equation for frequency of excited elec- 
trostatic waves: 



co 2 + ivco - k 2 Vj; -/3(kV se ) 2 coJ^ 



J 2 (b) 
co - joO 



= 0, (11) 



where: 

V = V in + Vei, W = joQ 



V se -- 

For Vi„ an d v e/ me expressions given in Braginskii (1963) are 
used. 

Because of the nonlinear character of the investigated pro- 
cesses, the averaging by casual phase the terms of hydrody- 
namic equations, which depend on *P do not vanish. 

In Eq. 9 we see that resonance occurs when co » joCl. It is 
convenient to represent co in the following way: 



CO-COQ + id. 



(12) 



here coq = joQ 

after substitution of Eq. 12 into Eq. 1 1 we obtain expression 
for frequency of the excited wave and their increment: 



coo = + ^ + i 2 ^,. + /; o W(W If ) 2 , (13a) 

S = - V -± ^+0.5J\(b)l3(kV se ) 2 , (13ft) 

We could solve equations 13a and 13b numerically. Fig. 1 
shows values of the function f w (k, jo), where f u (k, jo) = co if the 



Eqs.l3a-b have solution and /^(k, jo) — if Eqs.l3a-b do not 
have solutions. Fig. 2 represents values of the function f$(k, jo), 
where f s (k, jo) = Sif the Eqs.l3a-b have solution and f s (k, jo) = 
if Eqs.l3a-b do not have solutions. When solving equations the 
parameters, given in paper by Vernazza et al. (1981), are used. 

In our model the only source of excitation of electrostatic 
waves is the oscillation of neutrals, in the case of absence 
of acoustic waves i.e. V = 0, in the fourth term of Eq. 1 1 
jj o (b) _ Jo(0) 



, therefore we have: 



co 2 + iv ei co - k 2 V 2 Ti -p(kV se ) 2 = 0, 
the solution of Eq. 14 is: 



= -^± X \'4+k 2 Vl i+ P(kVse) 2 . 



CO 



(14) 



(15) 



Eq. 15 represent the damping wave with decrement and 

dispersion relation coo = ± y]-f + + P(kV se ) 2 , this result 
shows the fact that in the absence of acoustic waves the electro- 
static waves could not be excited and where the motion of the 
neutrals does not play a role the damping of propagating elec- 
trostatic waves will occur. 

In the regions of the solar atmosphere where the plasma ion- 
ization rate is high the motion of the neutrals does not influ- 
ence the dynamics of the charged particles. Due to collisions 
between the charged particles the damping of the electrostatic 
wave will occur and the energy of electrostatic waves could be 
transformed into the thermal energy of the solar atmosphere. The 
Fig. 3 shows dependence of the collision frequencies v ei and v en 
on the altitude in the solar corona, we see that at high altitudes v e , 
(which is responsible for wave damping) is higher than v en . For 
the Fig. 3 the solar atmosphere model given by Vernazza et al. 
(1981) is used. The dependence of the characteristic frequencies 
on the altitude in the solar atmosphere is given in Krasnoselskikh 
et al. (2010), where the model of the sun described in Fontenla 
et al. (2009) is used. 



3. Conclusions 

We investigated possibility of transformation of neutral gas mo- 
tion energy into the electromagnetic energy, in the solar atmo- 
sphere. We described pumping of the acoustic wave energy into 
the energy of electrostatic waves in the regions of the solar atmo- 
sphere where concentration of the neutrals is much higher than 
concentration of the charged particles. 

With increasing of the height in the solar atmosphere the ion- 
ization rate of the plasma increases i.e the role of the motion 
of the neutrals in the process of wave generation decreases and 
waves damp due tot he collision of the charged particles. The en- 
ergy of the damping waves could be transferred in to the thermal 
energy of the solar atmosphere plasma, at higher altitudes the 
temperature of the solar atmosphere plasma increases (Vernazza 
etal. 1981). 

We used two fluid approach and the parametric instability 
technique is used. We have determined that for particular param- 
eters anergy of acoustic waves of certain frequencies and certain 
wavelength could be pumped into the energy of the electrostatic 
waves. Figs 1 and 2 what kind of acoustic waves could contribute 
into the growth of the electrostatic waves. 

In this paper the problem was investigated using two fluid 
MHD approximation, we find it interesting to study transforma- 
tion of neutral gas energy into the electromagnetic energy using 
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kinetic approach, these approach could lead to more precise re- 
sults. 
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Fig. 1. Dependence of "filtered", excited waves frequency on wave number and jo- Following particular parameters are chosen: 
n 0n = 10 21 rrT 3 , T = 401QK. 




Fig. 2. Dependence of excited waves "filtered" increment on wave number and jo- Following particular parameters are chosen: 
n 0n = 10 21 trT 3 , T = 4010K. 
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h (km) 

Fig. 3. Dependence of the collision frequencies v, f -dashed line and v en -solid line, on the altitude in the solar atmosphere. 
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